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We compare a recent lattice determination of the nucleon distribution amplitudes with other 
approaches and models. We study the nucleon distribution amplitudes up to twist 6 in next- 
to leading conformal spin and we also investigate conformal d-wave contributions to the leading 
twist distribution amplitude. With the help of light-cone sum rules one can relate the distribution 
amplitudes to the form factors of the nucleon or the A*' ^ A transition at intermediate values of the 
momentum transfer. We compare our results with experimental data in the range 1 GeV'^ < < 10 

. - . GeV'^. Keeping in mind that we are working only in LO QCD and NLO-QCD corrections might be 

sizeable we already obtain a surprisingly good agreement for the nucleon form factors G"/, G^j, G\ 

^— . and and for the A'^ — > A transition form factor ratios Rem and Rsm- 

(N ■ 
u '. 

Oh. 

<t^ ; I. INTRODUCTION 

^: 

. The nucleon distribution amplitudes represent the universal non-perturbativc input to numerous exclusive reactions, 
' see, e.g., [l[ for an early review. Taking corrections up to twist-6 [2| into account we compare different non-perturbative 
(— I , methods to determine the nucleon distribution amplitudes, in particular lattice simulations [1, QCD sum rule 

O i' estimates 0, 0] and a phenomenological model Q . For asymptotically large values of the momentum transfer the 
' ] form factors can be expressed as a convolution of two leading-twist distribution amplitudes with a hard - pcrturbatively 
calculable - scattering kernel 0j H, B [l^j Hll, [H, [Hi • This approach (pQCD) is formally proven in the oo 
(-H ■ limit, and currently there is the consensus that pQCD is not valid at experimentally accessible values of the momentum 
transfer. In light-cone sum rules [13, [3 were worked out which relate the nucleon distribution amplitudes to the 
experimentally accessible form factors of the nucleon at intermediate momentum transfer. Form factors arc interesting 
quantities per se, since they encode information about the structure of the investigated baryon. This interest riscd a 

■ lot in recent years, in particular because new data from JLAB [TqI. [20l [2ll. 2^ for the well-known electromagnetic form 
[ factors of the nucleon contradict common textbook- wisdom. See, e.g., [23 1 for a review and further references. To our 

■ knowledge light-cone sum rules are the only theoretical approach to determine form factors at intermediate momentum 
' transfer that incorporate consistently the purely perturbative approach (pQCD). This was explicitly shown in the 
, case of the pion form factor [23| . If one calculates the light-cone sum rules for the pion form factor to leading order 

■ and next-to-lcading order in QCD one can show that the a^-corrcctions include the pQCD result in the oo 
0^ , limit. In the case of baryon form factors the pQCD result is expected to be included in the 0{a'^) corrections to the 

■ light-cone sum rule calculation. Currently only leading-order sum rules for the baryon form factors are known and a 
J> ' part of the NLO QCD corrections to the nucleon form factors. 

" '~j ■ The paper is organized as follows. In section 2 we introduce the concept of distribution amplitudes, in section 3 we 
rS ' collect QCD sum rule predictions for the nucleon distribution amplitudes and in section 4 we shortly explain the lattice 
^ . determination of the moments of the nucleon distribution amplitudes. All these approaches, including the numerical 
■ ■ ■ ' results, are discussed in section 5. The light-cone sum rule formalism is introduced in section 6 where we also give a 
short overview over the current literature on light-cone sum rules for baryonic form factors. In the next three sections 
we compare light-cone sum rule predictions with different models of the nucleon distribution amplitude for the form 
factors of the nucleon and for the — > A transition. In section 7 we use the nucleon distribution amplitudes including 
next-to- leading conformal spin contributions to determine the form factors, in section 8 we make use of some relations 
between twist-4 and twist-3 parameters and in section 9 we investigate the effect of the d-wave contributions to the 
leading twist distribution amplitude. We conclude and summarize our results in section 11. 

In the appendix we give for the first time the full expression for all nucleon distribution amplitudes up to twist-6 
including also the d-wave contribution for the leading twist distribution amplitude. 



II. THE NUCLEON DISTRIBUTION AMPLITUDES 



The distribution amplitudes comprise the infrared behaviour in exclusive processes involving large momentum 
transfer. They remove the infrared divergences in the perturbative diagrams encoding the nonpcrturbativc content 
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of the process and are defined in terms of the Bethe-Salpeter wave function 

*Bs(a;,fc±) - {0\T[q{xi,ki,^)q{x2,k2,±)q{x3,k3.±)]\P) (1) 

with Xi being the longitudinal momentum fraction carried by the quark i, ki^± its transverse momentum and \P) the 
nucleon state with momentum P [P^ ~ ^^n)- "^^^ distribution amplitudes are then obtained by integrating out the 
transverse momenta, 

^x^,^l) ^ Z{fi) d'h^^^Bsix^k^), (2) 

where Z results from the renormalisation of the quark field operators. In coordinate space the nucleon distribution 
am^itudes are derived from the following non-local nucleon-to- vacuum matrix element (here we follow the definitions 

{0\eijkul^{aix)[aix, aox]i' ^iU^p {a2x)[a2X, aox]j' ^jd'^ {a3x)[a3X, aox]k\k\P) ; (3) 

u and d are quark field operators, a, (3 and 7 are Dirac indices, while i,j and k are color indices; x is an arbitrary 
light-like vector, x'^ = 0, while the Ui are real numbers that fulfill ai +02 +03 = 1. The gauge-factors [x, y] are defined 
as 



[x,y] = Pexp 



1 

^<7 J dt{x-y)^A^{tx + {l^t)y) 




(4) 



where path ordering V is implied. They render the matrix element in Eq. ([3]) gauge invariant. In the following 
formulas we omit the gauge factors in order to simplify the notation. 

The leading-twist contribution to the nucleon distribution amplitudes has been determined long time ago including 
terms of next-to-next-to leading conformal spin, see, e.g., [l| for an early review. We will compare different determina- 
tions of the arising non-perturbative parameters in section 5. Currently the nucleon distribution amplitudes have been 
expanded up to contributions of twist 6 in Q and the corresponding non-perturbative parameters were estimated in 
0,11 with QCD sum rules and in Q from a phenomenological model. Some of these parameters were also calculated 
on the lattice d, [3j S H^, HE] . So-called x^-corrections (corresponding to deviations from the light-like separations of 
the quark fields in Eq. ([3])) to the leading twist distribution amplitudes were determined in [a. El, [13, Eal > they are 
formally of twist 5. 

Using the symmetry properties of the quark fields the matrix element in ([3]) can be expanded in twist as 

m^^Mal^yp{a2x)d';{a3x)\p) = ^(r3)L0(r4); , (5) 

i 

where r3/4 are certain Dirac structures and the Fi are distribution amplitudes, which can be expanded into eigenstates 
of conformal symmetry. This results in terms containing local operators. These local operators are associated with 
the moments of the distribution amplitudes, which are defined as 

1 

pn,n,n, ^ J_ f ^n, ^n, ^n, (^^ ^ , X3) . (6) 



F 



i.N 







Here Fi{xi, X2, X3) stands for a distribution amplitude and F^ pf for its normalisation constant, which is chosen such 
that F""" = 1. The integration measure is defined as 

Vx = dxi dx2 dx3 (5(1 — xi — X2 — X3). (7) 

Thus momentum conservation implies for the moments of the distribution amplitudes the relation 

Further details on distributions amplitudes (with complete expressions and definitions up to twist 6) are summarized 
in the appendices. 

For the nucleon distribution amplitudes isospin symmetry and the presence of two quarks of the same type implies 
that the number of independent distribution amplitudes is reduced compared to the general case. In particular. 
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the leading-twist nucleon distribution amplitudes can be expressed in terms of only one independent distribution 
amplitude which is usually taken as 



ip{xi,X2,X3) = Vi(a;i,X2,a;3) - Ai{xi,X2,X3) 



(9) 



and is equal to ^3{xi,X2,X3) in the notation of 0|- The distribution amplitudes Ai and Vi are defined in the 
appendices. At leading twist the nucleon distribution amplitude (p{xi) corresponds to the following form of the proton 
state [ii,!!!: 



P,T)= / Vx-^t^\uHxi)[u^ix2)d^ix3)-dHx2Wix3)]) 

Iq ^/9bXlX2X3 



(10) 



The first moments of ip{xi) can be interpreted as the momentum fractions carried by the quarks. 
The leading-twist distribution amplitude depends at leading conformal spin on one non-perturbative parameter, 
the normalization constant /jv, while for twist four we have two additional constants Ai and A2. In our approach 
no new parameters appear in leading conformal spin up to twist six. At next-to-leading conformal spin only two 
non-perturbative parameters Vi = (/s^"^ and A" = v?^*"^ — (^"^'^ arise in the case of leading twist and at next-to- 
leading twist we have three non-perturbative parameters, ff, /" and /2, for details see 0, 01. For the leading-twist 
distribution amplitude (p{xi) we have also determined the next-to-next-to leading conformal spin contributions which 
can be completely parametrized, e.g., by the moments v?^*^^, and The local matrix elements defining the 

non-perturbative parameters up to next-to leading conformal spin are (see for the corrected formulas from 



(0|£'^'= [u'C^u^ (0) b5^d%{0) 
(0|e*^'^- [u'Cj^u^ (0) b^rd^iO) 
(0|e'^"'= [u'Ca,,,uq (0) bs^^'d'^li (0) 
(0|£*^''= [u'C^u^] (0) [-f5^iz-Dd'']s{0) 
(0|£'-'''= [u'C^j5iz-Du^] (0) [^d%{0} 
(0| [u'Cj^u^] {0)[^^5Y'^z■Dd%{0) 



-ijk 



u'C-/^,-,5^z-D uA (0)[^7^d'^lA-(0) 



P) = Xim^NsiP) , 
P) = MmnNsiP) , 
P)^ fNVt{P-zf^N,{P), 

P)^-fNA1{P-zf^Ns{P), 
P) = \ift{P-z)M ^Ns{P), 
P) ^X2f2{P-z)M^Ns{P), 
P) ^\i.f^{P-z)M ^Ns{P), 



(11) 
(12) 

(13) 

(14) 

(15) 

(16) 
(17) 
(18) 



with the nucleon spinor Ns{P), the nucleon mass tojv, an arbitrary light-like vector z"^ with = and D — D- 
All derivatives act only on the quark fields and not on any explicit factor z. The second moments of the nucleon 
distribution amplitudes are related to the following local operators: 



[u'C^u^ (0) [75;^ iiz-D)'d%mP) = fN^'°'{P-zf^NsiP) , 



{{iz-DYu')C^u^ (0) [75^dl5(0)|P) 



(0|e*^''= \{{iz-Dfu')Ct^,v^] (0) [^d'=]5(0)|P) =/jv^'°°(P-z)V^5(P), 



.ijk 



{iz-Du')C^u^ (0) [-iz^iz-Dd%{Q)\P) 
{iz-Du')C^j5uA (0) iz-Dd'']s{0)\P) 



fNip'"\P-zr^Ns{P). 



(19) 
(20) 

(21) 



The parameters used in this work with their twist and conformal spin are summarised in the table below 





Leading twist 


Higher twist 


Leading conformal spin 
Next-to-leading conformal spin 
Next-to-next-to-leading conformal spin 


In 
A'l, Vf 

(plOl, v?200^ ^002 


Ai, A2 

J'u £d fd 
Jl J Jli J2 



As in the meson case these parameters can be estimated with QCD sum rules 13l| (see, e.g., [H, [H, [SJ] for some 
recent work in the meson case) or with lattice simulations (see, e.g., [ssl. [36j for lattice works considering the same 
mesons). 
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III. QCD SUM RULE DETERMINATION OF THE NUCLEON DISTRIBUTION AMPLITUDES 



The leading-twist distribution amplitude was investigated with QCD sum rules up to the second moments in 
[2^ . \3T\ and up to the third moments in [30| including perturbative contributions and terms proportional to the gluon 
condensate and to the four-quark condensate (several errors in [2^ were corrected in [sot). 

The next-to-leading twist normalization constants Ai and A2 describe the coupling to the proton of two independent 
proton interpolating fields used in QCD sum rules, Ai is the coupling of the so-called loffe current [3^, while A2 is the 
coupling of the interpolating nuclcon field that was advocated in [sl]. In 38 1 and [3§| first QCD sum rule estimates 
for Ai_2 were presented. Higher dimensional condensates were included in 40|. Unfortunately these pioneering works 
contain several misprints, for a review with the correct expressions see, e.g., (4ll. l42j. a^-corrections were calculated 
by Jamin in [i^. They turned out to be very large (« -1-50% for jA^I, corresponding to « 4-25% for |Ai|), but we 
will not take them into account, since we also do not have a^-corrections for the light-cone sum rules, connecting 
the distribution amplitudes with the nuclcon form factors. In Q also contributions of non-planar diagrams to the 
dimension 8 condensates were included. Putting all this together (for the first time) the QCD sum rule expression for 
Ai reads 



M%E^ 



Ml 
So 
Ml 



L- 



4 m5_ 
3 Ml 



where Mb is the Borel parameter, sq is the continuum threshold and 



Er,iso/M^) = 1- 



n-1 

E 

k=0 



A/2 



L = 



^s{Ml) 



(2^)2(:^£g'2 



0.55 GeV^ , 
~ 0.47 GeV^ 



{qgGq) 



0.65 GeV^ 



(22) 

(23) 

(24) 

(25) 
(26) 

(27) 

The 



We have neglected in Eq. (|22|) the small as-corrections to the four-quark contribution proportional to 
corresponding formula for A2 can be found, e.g^in @. QCD sum rule estimates for the defined in Eqs 
PS)) were first presented in 0] and updated in Q. The parameter set which is obtained by QCD sum rules will be 
called sum-rule estimate in the following, we use the numerical values from [sot for the moments of the leading-twist 
distribution amplitude and the values from @ for /^r, Ai, A2, A\, V-f, ff, /" and /2 . 
In our analysis wc use two related parameter sets which are based on the QCD sum rule determination: 

• Demanding that all higher conformal contributions vanish, fixes Vf , , f i and /2 , while the values for /^r, 
Ai, A2 are taken from the QCD sum rule estimates or from the lattice calculation. This parameter set will be 
called asymptotic. In the case of the leading twist, one would be left with the asymptotic distribution amplitude 
ip{xi,Q^ 00) ~ 'fiasyixi) = 12QxiX2X^fN ■ The corresponding expressions for the higher twist distribution 
amplitudes can be found in 0]. 

• With the help of light-cone sum rules [1, [l^, [13] one can express the nucleon form factors in terms of the eight 
non-perturbativc parameters /at, Ai, A2, A", V^^ ff, /" and /2 (inlcuding twist-6 corrections and expanding 
the distribution amplitudes up to the next-to-leading conformal spin). Choosing values for these parameters in 
between the asymptotic and the sum-rule values, we got an astonishingly good agreement with the experimental 
numbers, see |6|. This procedure is obviously rather ad- hoc and has to be replaced by a real fit after a^- 
corrections to the light-cone sum rules have been calculated. The paramter set obtained in @ will be called 
BLW. 



IV. LATTICE DETERMINATION OF THE NUCLEON DISTRIBUTION AMPLITUDES 



Lattice QCD offers the possibility to perform non-perturbative computations in QCD without additional model 
assumptions. For example, one can evaluate hadron masses and matrix elements of local operators between hadron 
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states. In particular, the non-perturbative parameters /tv, ■ ■ ■ introduced above can be extracted from Monte Carlo 
simulations on the lattice as advocated in [45| . 

Recently, the QCDSF collaboration has performed such a calculation It is based on gauge field configurations 
generated with two dynamical flavours of quarks. For the gauge field the standard Wilson action was used, while the 
lattice action for the quarks was the so-called non-perturbatively 0{a) improved Wilson fcrmion action, also known 
as the clover fermion action. Although lattice artefacts seem to be small, a reliable continuum extrapolation could not 
be attempted, and we utilize here the data obtained on the finest lattice corresponding to a gauge coupling parameter 
P = 5.40. Setting the scale via a Sommer parameter of ro = 0.467 fm the lattice spacing turns out to be a w 0.067 fm. 
On the lattice, matrix elements between the vacuum and a nucleon such as those needed here are computed from 
two-point correlation functions of the local operator Oa (x) under study and a suitable interpolating field Afa (x) for 
the nucleon. Asymptotically this two-point function decays exponentially with the distance between the operators 
since the lattice calculations are performed in Euclidean space. Projecting onto definite momentum one finds for 
sufficiently large (Euclidean) times t: 

^ ^ e-'^P-^e^P-y{0^ix, tWM 0)) = ^^^^^o {e{P - ^P ■ j + m^) e'^^^ )* . (28) 

Here Vg denotes the spatial volume of the lattice and the matrix elements of Oa {x) and Afa {x) have been represented 
as 

(O|0„(O)|P) =MoiV„(P), (29) 
(P|A/;(0)|0) = V^7V„(P). (30) 

As the local operators Oa{x) used in the simulations are linear combinations of the operators appearing in pT|) - (PT|) . 
the constants Mq are directly related to moments of the distribution amplitudes. 

The operators Oa{x) need to be renormalized. In Ref. [3| a non-perturbative renormalization procedure has been 
chosen. As the space-time symmetry on the lattice is reduced to the finite (spinorial) hypercubic group, the mixing 
pattern of our three-quark operators is more complicated than in the continuum and the choice of the operators 
becomes an important issue. Guided by the group-theoretical classification of three-quark operators given in [4^ the 
problematic mixing with lower-dimensional operators could however be completely avoided. Moreover, the freedom 
in the choice of the operators has ben exploited in order to reduce the statistical uncertainties of the results. 
Primarily, the combination of moments 

3 3 
has been evaluated, from which the combination (p"i"2n3 usually used in sum rule calculations is readily obtained by 

In the following sections we shall compare these lattice results with results obtained from other approaches and see 
what the lattice numbers imply for the nucleon form factors. 

V. COMPARISON OF DIFFERENT METHODS TO DETERMINE THE DISTRIBUTION 

AMPLITUDES 

In Table U we compare different estimates for the moments of the leading- twist distribution amplitude at 1 GeV^. It 
turns out that the BLW model, the BK model and the lattice evaluation give almost indentical results, which are close 
to the asymptotic values, while the QCD sum rule estimates seem to overestimate the deviation from the asymptotic 
form, although the deviation goes in the right direction. The BLW model was inspired by this experience: One starts 
with the asymptotic form and goes then in the direction of the QCD sum rule estimate, but only for a fraction of the 
whole difference. Choosing this fraction to be 1/3 one gets an astonishingly good agreement between light-cone sum 
rule predictions for the nucleon form factors a nd e xperiment, see In the same spirit one can make a BLW model 
for the second moments, also given in table IlldiSi . These values are again very close to the lattice values. The BK 



model [47[ was also inspired by experiment, in particular the decay J/"^ — *■ NN, the Feynman contribution to the 
nucleon form factor and the valence quark distribution function. In Table HIl we compare different estimates for Jn, 
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KS 


BK 
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1 
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5 

28 
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1 
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0.089 


0.12 


1 
8 


« 0.13* 


0.132 


1311 n 13'l('382'l 


^^^^ 


2 

21 


« 0.095 


-0.030(30) 


0.027 


0.02 


1 

12 


« 0.08* 


0.071 


0.0613(89)(319) 




2 

21 


« 0.095 


0.102(12) 


0.113 


0.10 


17 
168 


« 0.10* 


0.097 


0.1091(41)(152) 




2 

21 


« 0.095 


0.090(10) 


0.097 


0.10 


8 

21 


« 0.10* 


0.093 


0.1092(67)(219) 



TABLE I: Comparison of diflerent estimates for the moments of the leading-twist distribution amplitude renormalized at 1 
GeV'^. We show the asymptotic values (Asy) and the QCD sum rule estimates from [30| . Inspired by the QCD sum rule 
calculation two models for the leading-twist distribution amplitude were suggested, the COZ mode l KjOll and the KS model [33 . 
Using also some experimental input two phenomenological models were introduced, the BK model [431 ^-he BLW model (a]. 
Finally we show the lattice values from [3, 0, 0] . The first error is statistical, the second error represents the uncertainty due to 
the chiral extrapolation and renormalization. For the BK model no contributions from next-to-next-to-leading conformal spin 
were taken into account, thus the second moments denoted by the * do not contain any additional information and are fully 
determined by the first moments. 





As. 


QCD-SR 


BK 


BLW 


LAT 


/at • 10^[GeV^] 
Ai • 103[GeV^] 
A2 • 103[GeV^] 


5.0(5) 
-27(9) 
54(19) 


5.0(5) 
-27(9) 
54(19) 


6.64 


5.0(5) 

-27(9) 

54(19) 


3.234(63)(86) 
-35.57(65)(136) 
70.02(128)(268) 


Vf 




i « 0.333 


0.38(15) 
0.23(3) 


j- « 0.071 
If « 0.31 


0.13 
0.30 


0.1013(81)(298) 
0.3015(32)(106) 


ft 

A" 

fi 


0.30 
0.10 

4. w 0.267 

1 5 


0.40(5) 
0.07(5) 
0.22(5) 




0.33 
0.09 
0.25 





TABLE II: Comparison of difi'erent estimates for /jv, Ai, A2, A", Vi, fi, /" and /| renormalized at 1 GeV^. The QCD sum 
rule estimates and the BLW values are taken from [g| and we also show the phenomenological model from [ij] (BK). For the 
asympotic and the BLW parameters the values for /jv, Ai and A2 coincide with the ones from the QCD sum rule estimates. 

Ai, A2, A", Vi, fi, fi and at 1 GeV^. The leading twist parameters A", V-f, are already fully contained in table 
H] via the relations 

A^ = 2(^1"° -f - 1 , (33) 
= (34) 

Let us stress, however, that the errors quoted in Tables |T] and have to be taken with caution. On the lattice 
side a continuum extrapolation could not be attempted, and hence the associated systematic error is not included. 
Moreover, the errors on A" have been calculated by error propagation, which might not be too reliable. For the sum 
rule estimates the radiative corrections are expected to be sizeable, but these are only known for Ai and A2. The 
central lattice value of /at is about 35% smaller than the QCD sum rule estimates, while the lattice results for |Ai| 
and IA2I are about 30% larger than the QCD sum rule estimates. For Ai and A2 the discrepancy is reduced strongly, 
if radiative corrections to the sum rule estimates are included, cf. Eq. ([^^ . for Jn - according to our knowledge - no 
radiative corrections have been calculated yet. 

The parameters Ai and A2 can also be extracted from the lattice calculation of the nucleon decay matrix elements 
(expressed in terms of the parameters a and (3) in [2^, [2^ . Using the relations 



Ai = a , A2 = /3 , 



(35) 
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we obtain from the results in [2^, [2^ 

Ai = -43.90 ± 4.7 ± 8.5 • lO^^GeV^ , A2 = 93.96 ± 10.2 ± 22.7 • lO^^GeV^ , 

at the renormalization scale 1 GeV^. In that case the deviation from the QCD sum rule values is even more pronounced. 
In the non-relativistic limit one gets 

2Ai + A2 = 0. 

The estimates presented in table lU fulfill this relation almost perfectly: 

0± 0.24, (36) 

0.008 ±0.013. (37) 
For the ratio /Ai the differences between the central lattice and QCD sum rule estimates are even more enhanced: 

0.185 ±0.064, (38) 

0.0909 ±0.0054 ±0.0095. (39) 

The QCD sum rule estimate is a factor of two larger than the lattice result. In the next section we will see that the 
electromagnetic form factors of the nucleon depend only on the ratio but not on the individual values of and Ai, 
if the so-called loffe interpolating field is used, while the ^ A transition depends on the individual values. 



2Ai + A2 



2Ai - A2 
2Ai - 



QCD-SR 
- A2 



2Ai - A2 



LAT 



In 
Ai 



QCD-SR 

In' 



^1 / LAT 



VI. LIGHT-CONE SUM RULES FOR FORM FACTORS 

Light-cone sum rules (LCSR) are an advancement of QCD sum rules [3l| for intermediate values of the momentum 
transfer Q^, i.e., 1 GeV^ < < 10 GeV^ in the case of nucleon form factors. They were introduced in (TtI. [isj. The 
starting point is a correlation function of the form 

T(P,g) = J d'xe-^'^^0\T{r,iO)j{x)}\NiP)) , (40) 

which describes the transition of a baryon B with momentum P — q to the nucleon N{P) via the current j. The 
baryon B is created by the interpolating three-quark field 77. If i? is a nucleon one can use, e.g., the loffe current (sl] 
for the proton 

VioSeix) = e'^'^ [uHx)iCj.) uHx)] dUx) ■ (41) 

A typical example for j is the electromagnetic current in the case of the electromagnetic form factors 

= euu{x)jf,u{x) + edd{x)jf,d{x) . (42) 



With the definitions in Eqs. ([41]) . (|42)) the correlation function in Eq. (|40|) describes the electromagnetic form factors 
of the nucleon, which can be measured, e.g., in elastic electron-proton scattering. 

The basic idea of the light-cone sum rule approach is to calculate the correlation function in Eq. both on the 
hadron level (expressed in terms of form factors) and on the quark level (expressed in terms of the nucleon distribution 
amplitudes). Equating both results and performing a Borel transformation to suppress higher mass states one can 
express the form factors in terms of the eight (taking only leading and next-to leading conformal spin into account) 
non-perturbative parameters of the nucleon distribution amplitudes, the Borel parameter Mb and the continuum 
threshold sq, for details see d, [l^. 

We studied the electromagnetic nucleon form factors with the Chcrnyak-Zhitnitsky interpolating field {rjcz) in [iF 



In |44| we found that r]cz yields large unphysical isospin violating effects, therefore we introduced a new isospin 
respecting CZ-like current to determine the electromagnetic form factors. In Q we also studied the loffe current for 
the nucleon and extended our studies from the electromagnetic form factors to axial form factors, pseudoscalar form 
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factors and the neutron to proton transition. It turned out that the lofFe current yields the most rehable results. 
Despite our "bad experience" -qcz was used to determine the scalar form factor of the nucleon and the axial and 
the pseudoscalar one in [4§| . The question of the ideal interpolating field can also be addressed more generally: One 
can write down the most general interpolating field - without derivatives - of the nucleon as a linear combination of 
two currents and then try to optimize the relative strength of these currents. This approach was used for the scalar 
form factor of the nucleon in (50l | , for the axial- vector form factors in j5l| and for the electromagnetic form factors 
in [5^ . Since in (s^ x^-corrections were not included and different Dirac projections to extract the sum rules were 
used, we cannot easily compare the result with 

The light-cone sum rule method can also be applied to other observables than the nucleon form factors. 
In the class of nucleon to resonance transitions the following processes were considered: The TV — > A transition was 
studied in this framework i n IssI (for a similar approach for = see, e.g., [13]), the axial part of the N A 
transition was calculated in [55|. Very recently the form factors of the N — > A^*(1535) transition were presented in 
[5^. In [53, m, [5^ pion-electroproduction was investigated. 

Also decays of baryons can be described with that formalism: A;, — > plv was discussed in I27I . The authors of [60l | 
considered Ac — > Alf and therefore determined a part of the A distribution amplitude. In [6ll | the transition E — > A'^ 
was investigated. Recently the rare decays Ab — > A7 and Ab — s- Al~^l~ were treated in [6^ with the same formalism. 
Electromagnetic form factors of E and A-baryons were estimated in [gl, [g^l ■ 

So far all mentioned LCSR calculations for the baryon form factors were done in leading order QCD. One expects 
sizeable radiative corrections of up to 30 %. In [6^ a first step in calculating the full C'(as)-corrections to the nucleon 
electromagnetic form factors was performed. The intrinsic final uncertainty of this approach is expected to be in the 
range of less than ±20%, if QCD corrections are included. Comparing the theoretical predictions with experimental 
numbers one must be careful to distinguish between quantities directly calculated like Fi and F2 and quantities like 
Ge = Fi — (5^/(4m^)F2 for which cancellations might ruin the predictive power. 

In the following we use the LO QCD light-cone sum rules of for the electromagnetic form factors of the nucleon 
and the LO QCD results of [H^ for the N A transition to compare the consequences for the form factors which 
the lattice results for the nucleon distribution amplitudes entail with those which result from different QCD sum rule 
estimates. Note, however, that the errors on the non-perturbative parameters of the nucleon distribution amplitudes 
will not be taken into account, because this would not make much sense due to the inherent uncertainty in the LO 
light-cone sum rules. 



VII. RESULTS FOR THE FORM FACTORS AT INTERMEDIATE MOMENTUM TRANSFER 

In this section we use light-cone sum rules to extract physical form factors from the nucleon distribution amplitudes, 
by taking into account conformal spin contributions up to the p-wave; d-wave effects will be discussed in section 9. 
We compare our theory results to the following experimental numbers. For the electromagnetic nucleon form factors 
we take data from: 

• The magnetic form factor of the proton normalized to the dipole form factor G\J{i1pGd) from [66l. [stI. [gsI. [gqI. 

[zO,[Zi|,[22,[2i[2i, with 

Gd{Q'')^- -2, Mp = 2.7928... (43) 

(^+o.7iGcvO 

The data of [H, [6^, [T^I are actually taken from the reanalysis in [tHI ■ 

• The ratio of the electric and magnetic form factors of the proton /ipG^^/G?, from Rosenbluth separation [H, 
[63, [ei, [ei, [iS, [p, [p, [li [li, [li, O, [li and from polarization transfer [li, [lO, [lU, [H . We would like to point 
out here that [73, uM claimed already in the seventies a steeper dependence of compared to for 
momentum transfers above 1 GeV^. Currently the Rosenbluth separation data for Ge are judged to be less 
reliable. 

• The ratio of the proton form factors F^ and i^f given as y/O^Ff /{{fip — 1)-Ff) in [6^. [btI. F79| . 

• The magnetic form factor of the neutron normalized to the dipole form factor: G^j/dinGo) from [tsI. [soI. [Sll. 
[8i,[8l[8l[ll,[li with /i„ = -1.913... 

• The electric form factor of the neutron normalized to the dipole form factor: G^/G^) from [tsI. [soI. [ssI . [stI. [ssI . 
S2t M, i91, .92„ [3 [9I [96,, ,97,, ,98,, M. JM] ■ The data are very well described by the so-called Galster fit[m| , 
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we show in our plots the update of the Galster-fit from Kelly [102l |: 

f-,n,Galster(r^2\ (1.70 ± 0.04) T 2n ^ (aa-. 

^^)=l + (3.30±0.32)r^^(^)' ^^'^ ^ = 4^^ ' ^^^^ 



For the axial form factors we compare our result to the dipole formula [10; 



1 2fi7 

Ga{Q') = ■ (45) 

1^ QL 



(1.014GeV)5 



For more details see [6|. 

Finally we use the following data for the iV ^ A transition: 

• The magnetic form factor normahzed to the dipole form factor GJ^/(3Gd) from [loi, [Toi, [Toi, [lo3, [IqI [Tol 

[iTo[fiTi|. 

• The ratio of the electric quadrupole to the magnetic form factor Rem from |llll . lll j |. 

• The ratio of the Coulomb quadrupole to the magnetic form factor Rsm from [llll . Ill2l | . 

For more details see [s^. Since we compare the data with the LCSR predictions, which are expected to work best 
in the region 1 GeV^ < < 10 GeV^, we are only interested in experiments where values of the form factors for 
momentum transfer above = \ GeV^ are available. 

For the theory prediction we will use six models (i.e. six determinations for the nonperturbative parameters /at, Ai, 
A2, A\, Vf, fi, fi and ) for the nucleon distribution amplitudes including s- and p-wave contributions: 

1. QCD sum rule estimates (dotted red lines), 

2. asymptotic form (dashed red lines), 

3. BLW model (solid red lines), 

4. lattice evaluation plus QCD sum rule estimate for /| (dotted blue lines), 

5. lattice evaluation plus asymptotic values for (dashed blue lines), 

6. lattice evaluation plus BLW estimate for /| (solid blue lines). 

Since ff, and have not been determined on the lattice, we have to use in the lattice parameter set QCD sum 
rule estimates, asymptotic values or the BLW model for For the nucleon form factors we use the LCSRs obtained 
in @ and for the TV — > A transition we use the LCSRs obtained in . 

To our accuracy, the sum rules for the nucleon form factors depend only on the five parameters /jv/^i, A^, Vf, /" 
and /f . Within the light -cone sum rule approach we determine the form factors Fi and F2 directly. The electric and 
the magnetic form factors Ge and Gm are linear combinations of Fi and F2: 



GMiQ^)^Fi{Q^)+F2iQ^), 

Ge{Q')^F,{Q^)-^F2{Q'). (46) 



N 



As discussed above, Eq. (pS)) shows that in G^; cancellations occur. Therefore our predictions for Ge are less reliable 
than those for Gm- 

The light-cone sum rule predictions for the form factors are shown in Figs. [T] - [H For G\j,G^,G^j and G^ the 
differences between the lattice determinations and the other approaches (asymptotic, QCD sum rule and BLW) are 
smaller than the expected overall uncertainties, i.e., the pairs of parameter sets (1) - (4), (2) - (5) and (3) - (6) yield 
almost identical results. Since /" and ff were chosen to be identical within these pairs, differences can only occur 
due to Vi, Ai and /at/Ai. The lattice values for Vf and Ai are very close to the BLW values, while /at/Ai is in the 
lattice determination about a factor of two smaller than the QCD sum rule estimate. In Ge, \/Q^^2'/((Mp ~ l)-^f) 
and G^ cancellations occur, so we expect much bigger theoretical uncertainties and also big differences between our 
data sets might be possible. 

The data for G^j (Fig.[T]) are very well described with the asymptotic and the BLW data sets, the differences between 
the parameter sets (2) and (5) and between (3) and (6) are negligible. The pure QCD sum rule estimates (set (1) 
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FIG. 1: LCSR results for the electromagnetic form factors (left: Gm/{^pGd) vs. Q^; right: fipGE /Gm vs. Q^) of the proton, 
obtained using the BLW model (red solid line), the asymptotic model (red dashed line) and the QCD sum rule model (red 
dotted line) of the nucleon distribution amplitudes. The corresponding results for the lattice values of the nucleon distribution 
amplitudes are given in blue. The red data points on the right picture are JLAB data, while the blue and the green ones are 
obtained via Rosenbluth separation. Currently the Rosenbluth separation data for Ge are judged to be less reliable. 
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FIG. 2: LCSR results for the electromagnetic form factors of the neutron (left: Gm /(fJ-nGn) vs. Q^; right: Ge/{Gd) vs. Q^), 
obtained using the BLW model (red solid line), the asymptotic model (red dashed line) and the QCD sum rule model (red 
dotted line) of the nucleon distribution amplitudes. The corresponding results for the lattice values of the nucleon distribution 
amplitudes are given in blue. The thin solid blue line represents the updated Galster fit. 
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FIG. 3: LCSR results (solid curves) for the axial form factor of the proton Ga normalized to Gd = flA/(l + Q^/tt^a)^ vs. 
Q^(left panel) and the tensor form factor Gt normalized to Ga vs. (right panel), obtained using the BLW model (red 
solid line), the asymptotic model (red dashed line) and the QCD sum rule model (red dotted line) of the nucleon distribution 
amplitudes. The corresponding results for the lattice values of the nucleon distribution amplitudes are given in blue. 
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FIG. 4: LCSR results (solid curves) for the ratio \/Q^F2 /{F[ * 1.79) obtained using the BLW model (solid line), the asymptotic 
model (dashed line) and the QCD sum rule model (dotted line) of the nucleon distribution amplitudes. The corresponding 
results for the lattice values of the nucleon distribution amplitudes are given in blue. Red symbols: experimental values obtained 
via Polarization transfer. Blue symbols: experimental values obtained via Rosenbluth separation. Currently the Rosenbluth 
separation data for Ge are judged to be less reliable. 



and (4)) are about a factor of two too small. In the case of G^j (Fig. [5]) one sees the same structure for the different 
models of the nucleon distribution amplitude as for G^^, but now all theory predictions arc shifted to lower values. 
For (Fig. [3]) we agree for values below 5 GeV^ very well with the dipolc behavior, if we use the asymptotic or 
the BLW parameters; for higher we predict a slightly steeper fall-off. Again the pure QCD sum rule estimates are 
considerably worse. 

An interesting test of our approach is whether the unphysical tensor form factor Gt (Fig. [3]) is consistent with zero. 
This holds for all parameters sets except the pure QCD sum rule estimates (set (1)). In our approach G^ is not 
exactly zero, because we treat the intial proton state differently from the final proton state: one is described by an 
interpolating nucleon field, the other by the nucleon distribution amplitude. 

Finally we have the ratios G|/Gd (Fig. HD , G^/G^^ (Fig. [T|), and y/Q^pP/Hfip - l)Ff) (Fig. H]), which are very 
sensitive to the explicit form of the nucleon distribution amplitudes due to cancellations in Ge- If we just look at 
G^ our result would be consistent with zero and therefore describes the data well. If we investigate G^/Gd, we 
blow up the large contributions. Now we have a "perfect" agreement between the pure QCD sum rule parameters 
and the data. Our data set (4) is almost identical to the updated Galster fit. The BLW model is consistent with 
zero, while the asymptotic distribution amplitude yields negative values. The difference between the lattice values 
for the distribution amplitudes and the data sets (1) - (3) is visible, but not dramatic. In the case of G^/G^, and 
a/Q^Fj /(-Ff * 1.79) we can make similar observations. The purely asymptotic values lie above (below) the data for 

^b/^m iVQ^^2 /i^i * l-'''9))i the BLW data set moves the results in the right direction, but not far enough. Our 
data set (1) is completely off, because it predicts a very small value for Ff. Now we also have big differences between 
the lattice values of the distribution amplitudes and the pure QCD sum rule values. 

Taking into account that the full C'(Q!s)-corrections to the LCSRs are not yet available, already a rough agreement 
of our approach with the data is a success. For G^j, G^, GJj and G^ we get unexpectedly "good" results if we use 
the BLW form or the asymptotic form of the nucleon distribution amplitude. The corresponding lattice values (set 
(5) and (6)) give similar results. For G'%/Gd, G^eI^^i^ ^nd ^/Q^Fi/{Ff * 1.79) the BLW model Hcs in the right baU 
park, but we would prefer a model for the distribution amplitudes which lies in the middle between the asymptotic 
and the pure QCD sum rule estimate (BLW is closer to the asymptotic value). 

In the transition form factors of 7*A^ — > A all eight non-perturbative parameters appear, see (53| . The results are 
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shown in Fig. [5l As expected, now the differences between the parameter set pairs (1) - (4), (2) - (5) and (3) - (6) are 
more pronounced. In the case of the theory curves generally tend to be more flat than the experimental data. 
The form factors obtained with the lattice values for the nucleon distribution amplitude lie considerably above the 
data sets (1) - (3). Above « 3 GeV^ the asymptotic distribution amplitude and the BLW distribution amplitude 
are close to the data. The fact that Rem is close to zero is reproduced very well with the BLW parameters (sets (3) 
and (6)) and the lattice plus asymptotic values (set(5)), while positive values are obtained with the purely asymptotic 
form (set (2)). One gets a negative result with the QCD sum rule determination of the nucleon distribution amplitude 
(sets (1) and (4)). In the case of Rsm the differences are not very pronounced, all values are close to zero. Altogether 
one has to conclude that while all approaches give the correct order of magnitude none gives a really convincing 
description of all data. However in view of the fact that the systematic uncertainities are even more pronounced than 
for the nucleon form factors more could not have been expected. 
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FIG. 5: j*N A transition form factors {\eit:G*M /{3Gd) vs. Q^, middIe:J?EA/ vs. Q^, right: Rsm vs. Q^) in the LCSR 
approach [s^l obtained using the BLW model (solid line), the asymptotic model (dashed line) and the QCD sum rule model 
(dotted line) of the nucleon distribution amplitudes. The corresponding form factors based on the lattice values for the nucleon 
distribution amplitudes are given in blue. 



VIII. REDUCING THE NUMBER OF INDEPENDENT PARAMETERS OF THE DISTRIBUTION 

AMPLITUDES 



In [ll3j | the following approximate relations between twist-4 and twist-3 parameters were derived: 

fd = 3 1 /jv 

^1 10 6 Ai ' 

r = 

•'^ 10 6 Ai ' 



(47) 



Using these relations, we can express the nucleon form factors in terms of only three independent parameters, namely 
Vf and A". For the comparison with the data we show now only two models for the remaining three parameters 
of the nucleon distribution amplitude: 

(a) Lattice determination of the distribution amplitude - blue curve. 

(b) BLW model - red curve. 

We obtain the following values for /f : 





(a) 


(b) 


asymptotic 


BLW 


QCD-SR 


Jl 


0.11 


0.13 


0.10 


0.09 


0.07 


ft 


0.31 


0.33 


0.30 


0.33 


0.40 



which are compared with our previous estimates for /{* and ff. In this approach ff lies between the asymptotic and 
the BLW value; /{* is also close to the asymptotic or the BLW value, but its deviation from the asymptotic value is in 
the "wrong" direction. In Fig. [5] we show the electromagnetic form factors of the proton, in Fig. [7]the electromagnetic 
form factors of the neutron, in Fig. [8] the axial and the tensor form factor of the proton and finally in Fig. [9] the ratio 
of the form factors F2 and Fi of the proton. In Fig. [TO] we show the three A^ ^ A transition form factors. In all cases 
we obtain results which are very close the BLW results of the previous section, so it seems that the nucleon form 
factors are very sensitive to the values of Vf and A", while the dependence on /at/Ai, and /" is less pronounced. 
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FIG. 6: LCSR results for the electromagnetic form factors (left: Gm / (fipGD) vs. Q^; right: fipCs /Gm vs. Q^) of the proton, 
obtained using the BLW model (red solid line) and the lattice prediction (blue solid line) for the nucleon distribution amplitudes. 
In both cases the parameters /f are determined from the twist-3 parameters, cf. Eq. (47). The red data points on the right 
picture are JLAB data, while the blue and the green ones are obtained via Rosenbluth separation. 
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FIG. 7: LCSR results for the electromagnetic form factors of the neutron (left: Gm/{ij,„Gd) vs. Q^; right: Ge/(Gd) vs. 
Q^), obtained using the BLW model (red solid line) and the lattice prediction (blue solid line) for the nucleon distribution 
amplitudes. In both cases the parameters /f are determined from the twist-3 parameters, cf. Eq. (47). The thin solid blue 
line represents the updated Galster fit. 
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FIG. 8: LCSR results (solid curves) for the axial form factor of the proton Ga normalized to Gd = pyi/(l + Q'^/iti'a)^ vs. 
Q^(left panel) and the tensor form factor Gt normalized to Ga vs. (right panel), obtained using the BLW model (red solid 
line) and the lattice prediction (blue solid line) for the nucleon distribution amplitudes. In both cases the parameters /f are 
determined from the twist-3 parameters, cf. Eq. (47). 
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FIG. 9: LCSR results (solid curves) for the ratio sfO^F^IiFl * 1.79) obtained using the BLW model (red solid line) and the 
lattice prediction (blue solid line) for the nucleon distribution amplitudes. In both cases the parameters /f are determined from 
the twist-3 parameters, cf. Eq. (47). Red symbols: experimental values obtained via Polarization transfer. Blue symbols: 
experimental values obtained via Rosenbluth separation. 
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FIG. 10: 7* A'' A transition form factors {\%i\,:G\i j (?>G d) vs. Q^, middle:_R_BM vs. Q'^, right: Rsm vs. Q'^) in the 
LCSR approach [s^] obtained using the BLW model (red solid line) and the lattice prediction (blue solid line) for the nucleon 
distribution amplitudes. In both cases the parameters /f are determined from the twist-3 parameters, cf. Eq. (47). 



IX. EFFECTS OF HIGHER CONFORMAL SPIN CONTRIBUTIONS 



In this section we include (in comparison to the previous sections) also contributions of the next-to-ncxt-to leading 
conformal spin to the leading-twist distribution amplitude. These terms have been determined on the lattice 0, |3| 
and with QCD sum rules [29, 30l . [37j . The explicit expressions for the leading-twist distribution amplitudes can be 
found in appendix C. The contributions of these higher moments to the electromagnetic form factors of the nucleon 
have already been estimated in (l6j . but only for the Chernyak-Zhitnitsky interpolating field. Here we work out the 
contributions of the second moments to the light-cone sum rules for nucleon form factors using the loffe current. We 
will use the form in Eq. (jCl[) for the leading-twist distribution amplitude and the following parameter sets: 



1. Asymptotic distribution amplitude (black lines). 

2. BLW plus second moments from QCD sum rules (dotted red lines). 



3. BLW plus second moments from the lattice (dashed red lines). 
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FIG. 11: LCSR results for the electromagnetic form factors (left: GM/ifJ-pGn) vs. Q^; right: ^pGe/Gm vs. Q^) of the 
proton, obtained using the asymptotic form (black solid line), the BLW model (red) with second moments from QCD sum 
rules (dotted), lattice (dashed) and BLW (solid) and the lattice determination (blue) with from QCD sum rules (dotted), 
the asymptotic model (dashed) and the BLW model (solid). The red data points on the right picture are JLAB data, while 
the blue and green ones are obtained via Rosenbluth separation. 




FIG. 12: LCSR results for the electromagnetic form factors of the neutron (left: Gm/{^„Gd) vs. Q^; right: Ge vs. Q^), 
obtained using the asymptotic form (black solid line), the BLW model (red) with second moments from QCD sum rules (dotted), 
lattice (dashed) and BLW (solid) and the lattice determination (blue) with fl^ from QCD sum rules (dotted), the asymptotic 
model (dashed) and the BLW model (solid). 



4. BLW plus second moments a la BLW (solid red lines). 

5. Lattice evaluation plus QCD sum rule estimates for f^. (dotted blue lines). 

6. Lattice evaluation plus asymptotic values for (dashed blue lines). 

7. Lattice evaluation plus BLW estimates for (solid blue lines). 

In Fig. [TT]we show the electromagnetic form factors of the proton, in Fig. [T^lthe electromagnetic form factors of the 
neutron, in Fig. [T3] the axial and the tensor form factor of the proton and finally in Fig. [T3] the ratio of the form 
factors F2 and Fi of the proton. In Fig. [15] we show the three — > A transition form factors. 

In almost all cases the second moments of the leading-twist distribution amplitude determined with QCD sum rules 
give huge corrections. We show these parameter sets in the plots, but we will not discuss them any further. 
The magnetic form factor of the proton is very well described by the BLW model with second moments a la 
BLW (set(4)) or from the lattice (set(3)) and the lattice values for the distribution amplitude with /Jf from BLW 
(set (7)) or with the asymptotic values for /Jf (set (6)). This is not unexpected, since the second moments a la BLW 
and from the lattice are quite similar in size. This observation ensures, however, that there is not an unexpected 
strong sensitivity of the LCSRs to the second moments. The theory predictions for the magnetic form factor of the 
neutron C}^ are again shifted to lower values. Apart from this fact, the predictions for the parameter sets (1), (3), 
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FIG. 13: LCSR results (solid curves) for the axial form factor of the proton Ga normalized to Gd = <?yi/(l + Q^/o^a)^ "^s. 
(3^(left panel) and the tensor form factor Gt normalized to Ga vs. (right panel), obtained using the asymptotic form (black 
solid line), the BLW model (red) with second moments from QCD sum rules (dotted), lattice (dashed) and BLW (solid) and 
the lattice determination (blue) with from QCD sum rules (dotted), the asymptotic model (dashed) and the BLW model 
(solid). 
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FIG. 14: LCSR results (solid curves) for the ratio \/ Q^F^ l(F^ * 1.79) obtained using the asymptotic form (black solid line), 
the BLW model (red) with second moments from QCD sum rules (dotted), lattice (dashed) and BLW (solid) and the lattice 
determination (blue) with from QCD sum rules (dotted), the asymptotic model (dashed) and the BLW model (solid). Red 
symbols: experimental values obtained via Polarization transfer. Blue symbols: experimental values obtained via Rosenbluth 
separation. 




FIG. 15: 7*7V ^ A transition form factors (left:GM/(3G_D) vs. Q^ middle:7?i5A/ vs. Q^, right: Rsm vs. Q^) in the LCSR 
approach [53| obtained using the asymptotic form (black solid line), the BLW model (red) with second moments from QCD sum 
rules (dotted), lattice (dashed) and BLW (solid) and the lattice determination (blue) with from QCD sum rules (dotted), 
the asymptotic model (dashed) and the BLW model (solid). 
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(4), (6), (7) lie relatively close together and they agree a little better with experiment, compared to the case where 
the d-wave contributions have been neglected. Also for and we get nice results, unless we use the QCD sum 
rule values of the second moments. As expected, in G^/G^, G^/Gd and F2 /Ff cancellations arise that lead to a 
strong dependence on the concrete form of the nucleon distribution amplitudes. 

In the case of the — > A transition the inclusion of d-wave corrections leads to strong enhancements in the prediction 
of G^j, while Rem and Rsm agree now better with experiment. 

X. CONCLUSION 

We have compared a new determination of the nucleon distribution amplitudes based on lattice QCD with different 
values available in the literature. The non-perturbative parameters of non-leading conformal spin from the lattice 
evaluation turned out to be close to the asymptotic form and very close to the BLW model. For the leading conformal 
spin parameters /at, Ai and A2 the deviation between lattice and QCD sum rules is about 30%, which is possibly 
due to neglected radiative corrections in the QCD sum rules estimates. Our models for the nucleon distribution 
amplitudes can be related to measurable form factors with light-cone sum rules. Despite the fact that the light-cone 
sum rules are only calculated to leading order in QCD and despite an intrinsic uncertainty of light-cone sum rules of 
about ±20% we get a very good description of G^^, G^j, G^ and G^ at intermediate momentum transfer. In G^, 
G^ and F2 / cancellations occur, which limit our predictive power. In general we found the following tendency: 
The asymptotic distribution amplitudes describe the data already amazingly well. Pure QCD sum rule estimates for 
the non-perturbative parameters overestimate the deviation from the asymptotic form, but the deviation goes in the 
right direction. The best results are obtained for the BLW values and the very similar lattice values. Including also 
d-wave contributions to the twist-3 distribution amplitude improves the description of G^, G^, G^ and G^ a little 
bit, but results also in bigger uncertainites in G^, G^ and i^^/Ff . In the case of the TV ^ A transition, we do not 
see the steep fall-off of G\,j, but the smallness of Rem and Rsm is very well reproduced. 

Further improvements on the theoretical side can be achieved by determining the NLO QCD corrections to light-cone 
sum rules, which connect the nucleon distribution amplitudes to the form factors. To match the NLO QCD accuracy 
also cts-corrections have then to be included in all QCD sum rule estimates of the moments of the nucleon distribution 
amplitudes. 
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APPENDIX A: NUCLEON DISTRIBUTION AMPLITUDES UP TO TWIST-6 

For completeness we give in this appendix the fuU expressions for the nucleon distribution amphtudes up to twist-6, 
details can be found in [3, Q . The general Lorentz decomposition of the matrix element defined in Eq. ^ reads 0] 

4 (0| e'J'^<(ai.T)u^^(a2x)d^(a3a;) \P) 



Vi 



7 



+V^ml {^C)^p (75^)^ + V^ml h^C)^^ {la^''' x.^^N)^ + Vem% {^C)^^ {^l^N)^ 
+ {r^ + ^!^r/^^ (7^75iV)-, + TamAT (x^P''za^,C)„^ (75^^)^ 



7 



7 



(Al) 



with 







Si 


= ^1, 


2(P-x)52 = 


5*1 — 5*2 , 








Vi 




2{P-x)V2 = 


--P2-P1. 








Vi 




2{P-x)V2 = 


Vi~V2~V3, 








2V3 




4(P-x) V4 = 


-2Vi + V3 + Va + 


2V5, 


4(P 








4(P-a:)'V6 


= -Vi + ^2 + V"3 4 










= A, 


2(P-x)A = 


= -Al + ^2 - A3 , 








2^ 


= A, 


4(P-x)y^4 = 


= -2Ai - ^3 - A4 - 


f 2yl5, 






-A, 


4(P-x)'A 


= yli - yl2 + A3 + 


A4 - A5 + Ag , 






Ti 




2{P-x)T2 = 


Ti + T2 - 2T3 , 












2{P-x)% = 


T1-T2- 2T7 , 




X5 = 


-Ti- 




h2r8, 


4{P-xf% 


= 2T2 - 2r3 - 2P4 


+ 2T5 + 2r7 + 2T8 


4(P 


x)r7 


= T7 




^{P-xf% 


= -Ti + T2+T5- 


Tg + 2T7 + 278 ■ 



(A2) 

The calligraphic notation is used for distribution amplitudes belonging to a simple Dirac structure, while the 
non-calligraphic functions denote distribution amplitudes of definite twist. Each distribution amplitude F = 
Vi, Ai,Ti, Si, Pi can be represented as 



F{ai,a2,a3,{P-x)) = e'^^^'^) ^'"'P(x,) , 



(A3) 



where the functions F(xi) depend on the dimensionless variables Xi, Q < Xi < l,X]i-'^'' ~ ^ which correspond to the 
longitudinal momentum fractions carried by the quarks inside the nucleon. 
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APPENDIX B: EXPANSION OF THE NUCLEON DISTRIBUTION AMPLITUDES UP TO NEXT-TO 

LEADING CONFORMAL SPIN 

In [2| the distribution amplitudes were expanded up to next-to leading order in the conformal spin. The twist-3 
distribution amplitudes read 

V^{x,,^i) = 120x12:23:3 [03(m)+0J(A')(1~ 3x3)] , 

Ai{xi,^j) = 120X1X2X3(X2 - Xi)03 (/x) , 

Ti(x„^) = 120x1X2X3 ((/-J -(/>^)(At)(l- 3x3)] . (Bl) 

The twist-4 distribution amplitudes read 

V2{x.„y) = 24x1X2 [</)Jl(Ai) + 0j(M)(l -5x3)] , 

A2{Xi,y) = 24X1X2(X2 - Xi)04 (/x) , 

T2(x„Ai) = 24x1X2 (1-5x3)] , 

V3ix„^i) = 12x3 [i;li^i){l~x3) + ^J4i^I){^-x3- 10x1x2) 

+1p^{fl){xl + X2 - X3(l - X3))] , 

A3{x,,fi) = 12x3(x2-xi)[(V'2 + V'4)(Ai)+V'r(A')(l-2x3)] , 
T3(x„Ai) = 6x3[{4>l + ^JJ"^+eMil-X3) 



+ + V^l + e|)(Ai)(l -X3- 10x1x2) 



'^4 



+ (^4 -V'4 +C4 )(m) (3^1 +2^2 - 2:3(1 - X3))] , 
T7(x„m) = 6X3 [(<^^+V4-^2)(m)(1-2^3) 

+{4't + i^t - ^I)(aO(i - - 10x1x2) 
+(9^r - V'4r - 0(M)(a;? + 2;2 - 2^3(1 - 2:3))] , 

5i(x„m) = 6x3(x2-xi) [(0^ + v!^+e2 + </'t + ^-1 + e4+)(M) 

+ ((^4 -^4 +C)(m)(1-2x3)] , 

Pi(x„ m) = 6x3 (xi - X2) [(0^ + v!^ -i2 + ^t + y^t - (t)if^) 

+ (<^4 -V4"-0(M)(l-2a:3)] . (B2) 
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The twist-5 amplitudes are given by 

V4{x,,^l) = 3[^°(/i)(l-a;3)+V5+(Ai)(l-2;3-2(x?+x2)) 

+ ^'5'(m) (2a:;i2;2 - 2:3(1 - 2:3))] , 

3(2:2 - xi) [-V'5(m) + ^5^(Ai)(l - 2x3) + V^5"(m)2:3] 
3 



T4{x^,fi) = ^[(0^ + V'5°+e5)(M) (1-2:3) 



+ (05+ + ^+ + 4+) (/i)(l -X,- 2{xl + xD) 

+ (05^ + (m) (2X1X2 - 2:3(1 - 2:3))] 

n{x,,f^) = ^[(0^ + ^°-e5")(M)(i-2:3) 

+ (05^ + i't - ^t) (m)(1 - 2:3 - 2{xl + xD) 

+ ('?^5" - V's" - O (m) (22:1X2 - X3(l - X3))] 

V5{x,,^i) = 6x3 [0g(//) + 0^(Ai)(l-2x3)] , 
^5(xi,^) = 6x3(x2 - xi)(?:)jr(/i) , 
T5(xi,M) = 6x3 fe°(M) +4+(m)(1 - 2x3)] , 

^2(x„m) = ^(x2-Xi)[-((/)°+V5°+4°)(A^) 



+ 



*^+^'5''+e5'')(M)(l-2x3) 



P2(x„m) = -(x2-Xi)[-(-0°-V^° + e5°)(M) 

+ (-0^-v^^ + e^)(M)(i-2x3) 

+ (-05-+7/'5-+e5")(Ai)x3] , 

and the twist-6 contributions are given by 

Ve{x,,^i) = 2[4>",{^l)+^+{^l){l~3x3)], 
Ae{xi,^j,) = 2{x2 - xi)(j)Q , 

n{x,,^i) = 2Uo(M)-i(0+- 06-) (1-3x3) 



(B3) 



(B4) 



The coclBcients (f>f,tpf and (i stands for the twist) in the above expansions can be expressed in terms of the eight 
non-pcrturbative parameters /jv, Ai, A2, /2, ^1, Vi, defined in section 2. The corresponding relations read, for 

the leading conformal spin: 

0!l = 0g-/jv, 0^ = 0^ = ^(/w + Ai) , 
1 



^4 — ^5 

For the next-to-leading spin, for twist-3 



^4" = 4° = ^A2, = € 



6 



(/w-Ai) 



for twist-4: 



9a 



i[/Ar(3-10n^) + Ai(3-10/f)] 



(B5) 
(B6) 



[/Ar(l-2Ar)-Ai(l-2A^-4/r)], 
[/w(2 + 5A5'-5Vf ) - Ai(2-5/f-5/n] , 



[/at (2 - Ai" - 3^1^^) - Ai(2 - 7ft + /i")] 



16 



A2(4-15/2''), ^4 



16 



A2(4-15/2''), 



(B7) 
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for twist- 5: 









-Ai(l-4/{')], 


05" = 


-| [/a^(1 


- 2Ai") 


-Ai(/f-/r)], 


i^t = 




F2A5'- 


2T/i'^)- Ai(l-2/f-2/i" 




\ [fN{2 






it = 









(B8) 

OU 4: 

and for twist-6: 

-\ [fN{l - ^Vt) - Ai(l - 2ft)\ , (B9) 
06 = i [/a. (1 + 4A'i') + Ai (1 - 4/f - 2/r)] . 



Next wc summarize the expressions for the a;^-corrections to the leading twist distribution amplitudes Vi, Ai and 
Ti . These corrections have been determined in [3, [23, [1^ . 
For Vi we have 

1— 2:2 

f 2 
^"^(a;2) = / dx,V,^\xuX2, 1 - x, - X2) = ^ {InC^ + XiC^) , 





1—X3 



Vf('^(x3) = J dx,V,''{xul~x^-X3,xs) = ^{fNCf + XiCt) 



(BIO) 



with 



CJ = {1- X2f [113 + 495.T2 - 552x1 - 10A5'(1 - 3x2) + 2^1^^(113 - 951x2 + 828a;^)] , 
= -(l-X2)'[l3-20/i^ + 3x2 + 10/i"(l-3x2)], 
= -(l-xs) [1441 + 505x3- 3371^:3 + 3405x^- 1104x3 

- 24yi'^(207-3x3-368x^+412x3 - 138x;^)] - 12(73 - 220^1^') ln(x3), 

= -(1 - X3) [11 + 131x3 - 169x2 + 63x^ - 30/f (3 + IIX3 - 17x2 + 7x^)] 

-12(3-10/f)ln(x3). (Bll) 



In the case of one finds 

/'2 
dxi^f (Xi, X2, 1 - XI - X2) = ||(1 - X2)3 {InD]. + X,Dl) , 


1-X3 

Af'^''\x3) = f dxiAf (xi,l-a:i -X3,X3) = 0, (B12) 



1—X2 



with 



D^f = 11 + 45x2 -2A5'(113 - 951x2 + 828x2) + 10^^(1 -30x2), 
D"^ = 29 - 45x2 - 10/i"(7-9x2)-20/f(5-6x2). 



(B13) 
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Finally, for Ti one has 



1~X2 



with 



^1 {X2) 



^1 [xs) 



J dxiTf (zi, X2, 1 ~ XI - X2) = g {/nE^ + XiEl 





1-x 



{fNEj + MEi) 



587a;? 



IMxi) 



E't ^ - 



E". = 



(1 - X2) (3(439 + 71x2 - 621x^ 

4A5'(1 - X2)^(59 - 483x2 + 414x^) -4Vf (1301 - 619x2 - 769x^ + llGlxf - 414x^)) 

12(73- 220Vf)ln(x2) , 
(1 - X2)(5 - 211x2 + 281x^ - lllx^ 

+ 10(1 + 61x2 - 83x2 ^ 33a;3)/d _ 40(1 _ ^2)2(2 - 3x2)/nj - 12(3 - lOff) ln(x2) , 
17 + 92x3 + 12(^5" + V{^){3 - 23x3) , 
-7 + 20/f + 10/1". 



(B14) 



(B15) 



APPENDIX C: EXPANSION OF THE NUCLEON DISTRIBUTION AMPLITUDE OF TWIST-3 UP TO 

NEXT-TO-NEXT-TO LEADING CONFORMAL SPIN 

The expansion of the leading-twist distribution amplitude in a basis which is diagonal with respect to one-loop 
renormalization reads up to next-to-next-to-leading conformal spin 

2 r 8 
i^(xi,X2,X3,^) = 120xiX2X3/jv(a<o)-^^|1 + hini^J.o)ixi - 2X2 +X3)L^ 

+ /lll(/io)(xi-X3)L^ 

+ /i2o(/^o) [1 + 7(x2 - 2x1x3 - 2x^)] 

40 

+ /i2l(/^o) (1 - 4x2) (Xi - X3) L^f'o 

+ /i22(/^o) [3 - 9x2 + 8x^ - 12x1x3] } (CI) 



with 



/?o = 11 - . 



(C2) 
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The coefficients hij can be expressed in terms of tlic moments by 

hM^l{l-3v"'"{fi)) (C3) 

= -J(1-3K(m) + ^Aa^))) (C4) 

= -l{4>tif^)-4>3if^)) , (C5) 

= (C6) 

= ^{l + AU^^)-3V,'i^,)) (C7) 

= l{s^tip) + i3if^)) , (C8) 

/i2o(m) = y (hioif^) +4-7 (3^i"1(a.) + v^2"°(m) + ^""'(m))) , (C9) 

(m) = 126 (^200 (^) _ ^002 (^)) _ g;^^^ ^ (CIO) 

/i22(/^) = ^ i-hioip) -4 + 6 + 2^200(^) ^ 2/"2(^))) . (cii) 
o 

Of course, this fo rm of ip is not imiquely determined by the moments. The anomalous dimensions were obtained, e.g., 
in [TTi . [Tll[Tl6j . 

One can also write down the rcnormalization group equations for the moments (^"i"2n3 — y^"i»2n3 _ ^nin2n3 qIqj^q- 

{7 + hu>ipo)Lm + huML^) , (C12) 

=^ (7 - 2h,oMLm^ , (C13) 

^°''H^^) (7 + hioML^ - hnML^) , (C14) 



V>^°H^^) (12 + 3/iio(Mo)i^ - 2/i2o(/io)i^ - h22ML^) , (C15) 

=^ (36 + 9/iio(Mo)i* + 9/iii(/io)i^ 

14 40 32 \ 

+ h2o{^JLo)L^^o +/l2i(Aio)i'"o +3/l22(Mo)i'"'« J , (C16) 

</'°'"(m) =^ (36 + 9/iio(//o)i^ - 9/iii(/io)£^ 

14 40 32 \ 

+h2a{^lo)L■'^^> - h2i{jio)L'>t'» +3/i22(Mo)i"''« j • (C17) 

Next we determine Vi, ^1 and Ti from up to d-wave contributions. Including all anomalous dimensions we obtain 
(in the following we suppress the explicit rcnormalization scale dependence in the formulas) 

Vi{xi,X2,x^) = 120xiX2X3/ArL^ |l - ^(1 - 3x3)^^ + ^(1 - Sxa)^^ 

hor\ 14 
- [(-2 + 7(X1 +X2'- ^XiX2)] 



- ^ [-1 + 8x2 - 8x1 -X3~ &X2X3 + 4:xl] 



2 

^22 

2 



[6 - 21x1 + 20a;? - 21^2 + 24x1X2 + 20x3] L'^^o | , 



(CIS) 



2 r 8 20 

Ai{xi,X2,X3) = — 60xiX23:3(xi — X2)fNL^ I ShioL'^ + hiiL'^ 

( 14 40 32 \ 1 

7h2oL^^a +/i2iL^ +h22L^^an , (C19) 
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Ti(a;i,a;2,a;3) = 120xiX2X3fNL^ {l + /iio(l - 3.T3)iW 



/l20 [-1 + 14X1X2 - 7X3 + Uxl] 



3 + 12xiX2 + 9x3 - 8x^1 L 



32 'I 



(C20) 



In the light-cone sum rule determination of the nucleon form factors we need the distribution amphtudes 
renormahzation scale fi ~ hq, therefore we give also the simplified expressions {L = 1) in the following. 
Our expressions agree up to next-to- leading conformal spin with the correspondings ones of 0. In 
function Vi including the next-to-next-to-leading conformal spin was used with the following notation: 



at a certain 
the vector 



yi(xi,X2,X3) ^ 120xiX2X3/Ar 



Ai(xi,X2,X3) = 120X1X2X3(X2 - Xi)/jv 1^3 + ( 

1 ,~ 



1 + 4>tMi^ - 3a;3) + [3 - 21x3 + 28x2] 
+ 4,f[5{xl+xl)-3{l-X3f)] 

}■ ■ 



1 " 4x3)0f 



ri(xi,X2,X3,^) = 120xiX2X3/Ar <j 1 + ^((^3 - 4>t)i^ - ^^s) 



Td4 



Ids 



[1 - 14xiX2 + 7x3 - 14x3] 



1 -|- X1X2 — 8x3 + llx 



(C21) 
(C22) 



(C23) 



with 



If 



If 



— (/120 - h2i + 3/122) 
^(3 + 28/"2-2iyi'^) , 

^(-7/120 + 2/121+4/122) 
5 

_^ (3 + 5^- _ V,^ - 2((^""2 ^ 5^101 ^ 5^200)^ 

O 

\ (7/120 + /121 + /122) 



20 



t22 



= (3 + 7(^1' + V^) - 56(^"°2 _ 2^101 ^ ^200)^ 



If = 2h 



22 
63 
10 



(3 -f- 7{A-^ + V,^) - 8(2(^002 ^ ^101 ^ 2(^200)) 



(C24) 

(C25) 

(C26) 

(C27) 

(C28) 

(C29) 
(C30) 
(C31) 

(C32) 
(C33) 



Numerically one obtains for the COZ model [3C 



4>f{fi = IGeV) = 0.61, 
4>f{fi = IGcV) = 3.7. 



(C34) 
(C35) 



This agrees with the numbers quoted in [l^. Using the lattice calculation we get 



(Ijfifi = IGeV) = 0.51 
4>fifi = IGeV) = 0.71 



(C36) 
(C37) 



Here again the pure QCD sum rule calculation seems to overstimate the effects. 
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APPENDIX D: MODELS FOR THE LEADING-TWIST NUCLEON DISTRIBUTION AMPLITUDE 



111 this section we present concrete models for the leading-twist nucleon distribution amplitude including next-to- 
next-to-leading conformal spin at the renormalization scale 1 GeV. At twist-3 one independent distribution amplitude 
(f{xi,X2,X3, fi) arises, see, e.g., [3]: 

ip{xi,X2,X3,fi) = {Vi - Ai) {xi,X2,X3.,y) . (Dl) 

In this distribution amplitude was denoted by <f>3(a;i, a;2, a;3, /i). From Lp one easily gets Vi, Ai and Ti, see, e.g., 
i: 

Ti{xi,X2,x-i) = i [ip{xi,x:i,X2) + ip{x2,xz,xi)] , (D2) 

Vi{xi,X2,xz) = i Yp{xi,X2,xz) + 'p{x2,Xl,Xz)\ , (D3) 

Ai{xi,X2,X3) = - [ip[x2,Xl,X3) - ip[xi,X2,Xz)] . (D4) 

The asymptotic form - only the leading conformal spin contribution - of </'(a;i , j;2, xa, /i) reads 

ipAsy{xi,X2,Xz,tj) = l20xiX2X3.(i)l{ll), = f n) ■ (D5) 

Including next-to leading conformal spin one gets 0] 

ip{xi,X2.,Xz.,tJ.) ^ ipAsy{xi,X2.,Xz,^i) 1 -f 03 (^) (xi - X2 ) -|- (^j^ (^i) ( 1 - Sxs) , (D6) 

with 



h 1+ _ (t>t 



(D7) 



In the literature also second moments of the leading-twist distribution amplitude were determined with QCD sum 
rules (29l. [sol. [37j . With this information one can build models for (^(xi, 2:2, 2^3) at a certain renormalization scale /i, 
inclu ding next-to- next-to leading conformal spin. We will use the model from [sot Lp^^^ {xi,X2,X3) and the model 
from [331 {xi,X2tXz)'- 

ip^'^^ {xi,X2,Xz) = ipAsy{xi,X2,X3) 

[23.8142;? + 12.978x^ -f 6.174a;^ -I- 5.88x3 - 7.098] , (D8) 

Lp'^'^{xi,X2,X'i) = LpAsy{xi,X2,X3) 

[2O.I6X? -f 15.12x^ -f- 22.682;^ - 6.72x3 + 1.68(xi-X2)- 5.04] . (D9) 

Bolz and KroU derived a very simple model using some experimental constraints (47| . Their model for the leading- twist 
distribution amplitude reads 

Lp^^{xi,X2,X'i) = ]^(pAsy{xi,X2,X3){\ + Zxi) . (DIO) 

Based on the lattice calculations of we have obtained the model 



(/j-^-^^ (xi, X2, X3) = ipAsy{xi,X2,X3){ - 0.401 + 29.214xi - 44.542x2 + 7.664x: 



3 



12.561XX2 -I- 31.748x1X3 - 103.09x2X3 (Dll) 
41.880xi 4- 92.958x^ + 17.836x3) ■ 
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